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We present a theoretical approach which maps polymer blends onto mixtures of soft-colloidal 
particles. The analytical mesoscale pair distribution functions reproduce well data from united 
atom molecular dynamics simulations of polyolefin mixtures without fitting parameters. The theory 
exactly recovers the analytical expressions for density and concentration fluctuation structure factors 
■ of soft colloidal mixtures (liquid alloys) . 
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I. INTRODUCTION 

oo ■ 

Polymer blends are systems of fundamental scientific interest. Their structural and dynamical properties change as 
r- — i. a function of the proximity to thermodynamic conditions of phase separation (i.e., their spinodal curve). Blends have 
many practical applications since new materials can be produced with specific physical and chemical properties by 
simply mixing polymer melt components that have the desired characteristics. Finally, polymeric materials used in 
daily life are usually mixtures of polymers with different chain length and/or different local structure. Understanding 
how the mixing of polymer melts modifies their properties has been a relevant and longstanding scientific and techno- 
logical goal in polymer physics and engineering. [lj,|3,|3| From a theoretical perspective, relevant work in this direction 
has been developed by several groups in recent yearsTP, IE IE IE IE 

A great deal of information on the correlation between local (intra- and intermolecular) structure and global 
fluid properties has been obtained by computer simulations of polymer mixtures. [E El El 0, El El O ne of the 
challenges in simulating polymer blends is the large range of length- and timescales that need to be investigated. 
Blend properties strongly depend on the local bond lengthscale (short timescale), since mixtures of polymers with 
different local chemical structure can lead to dramatically different physical properties, including, for example, an 
£h ' opposite trend in demixing with temperature. On the other end, properties need to be investigated on the large 
lengthscale (long timescale) of concentration fluctuations, which diverges approaching the spinodal decomposition. In 
this regime, the box size of the system simulated becomes the upper limit in resolution. 

To expedite computer simulations, it is useful to renormalize the system by adopting a coarse-grained description of 
the liquid. The goal is to have a formalism that allows us to simulate macroscopic behavior while retaining information 
about the detailed atomistic chemical structure of polymer chains. The overall behavior of blend materials has to be 
related, through a well-defined and reversible procedure, to the local chemical structure of the system. As a first step, 
a convenient approach is to introduce a united atom (UA) description of the polymer chain: a commonly adopted 
method in simulations of polymer melts and blends. |l2l Il3l Il5l A more drastic renormalization at the lengthscale of 
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the polymer size was proposed by Dautenhahn and Hall|l6| and later on by Murat and Kremer[T^. who mapped a 
blend of polymer chains onto a liquid of colloidal particles interacting through phenomenological soft-core potentials. 
Recently, Hansen and coworkers have developed a rigorous numerical methodology to derive an effective potential for 



polymers in solution using liquid state theory. [Isl fx£ 

An analytical solution of the c.o.m. mean-force potential for a melt was derived by one of the authors. |2(J, |21j 
The potential is qualitatively in agreement with the physical behavior observed in simulation data. For example, 
when compared with the classical solution by Flory and Krigbaum, 22] it has the advantage of correctly predicting 
greater interpenetration between two chains with increasing degree of polymerization, chain stiffness, and/or liquid 
density. However, since the derivation does not include explicitly chain connectivity, the analytical solution reproduces 
the short-scale simulation behavior after fitting the value of the intermolecular potential at complete interpolymer 
overlap. |2l) More recently we derived an analytical expression for the potential of mean-force in polymer melts, which 
reproduces well, and with no adjustable parameters, data from UA Molecular Dynamics (MD) simulations of chains 
with different degree of polymerization and local chemical structure. [2^] This coarse-grained potential is an explicit 
function of atomistic chain parameters, such as degree of polymerization and local semiflexibility, thereby correctly 
bridging different lengthscales of interest. 

It is clearly advantageous to have an analytical form of the potential versus a numerical one, obtained form the 
inversion of computer simulation data. Since the mean-force potential is a free energy, it depends on the thermody- 
namic state of the system considered. A numerical calculation of the potential would require us to perform numerical 
simulations for each thermodynamic state of our system, defeating the purpose of adopting a coarse-grained model 
to reduce c omp utational time. In this spirit, Krakoviack et al. recently derived an analytical potential for polymers 
in solution. |24| 

In this paper we extend the melt theory to derive an analytical expression for the effective interaction potential 
acting between the centers of mass (c.o.m.) of two polymer chains in a blend. We start from a first-principles liquid 
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state description, thus limiting our theory to the miscible region of the phase diagram. However, it includes the 
buildup of concentration fluctuations as the system approaches its spinodal decomposition. It is known that the 
detailed nature of single-phase blend correlations (pair distribution functions) is sensitive to system-specific factors 
such as liquid density, temperature, blend composition, and differences in effective unit size and local architectural 
details of the components. Moreover, many blends of practical relevance are characterized by strong asymmetries 
in local chemical semiflexibility and architecture. These effects influence local entropic packing of molecules and 
can lead to phase separation upon heating in a blend characterized by a lower critical solution temperature. Our 
analytical expression describes the interaction betweeen c.o.m. of a pair of polymers in a blend as a function of all 
the characteristic parameters that govern blend structure and dynamics. 

The derived mesoscale pair distribution functions effectively map the polymeric liquid onto a fluid of interacting 
soft-colloidal particles, thereby correctly reproducing fluctuations in number density and concentration, i.e. struc- 
ture fact ors, and recovering the well-known equations for the compressibility and the dilation factor of colloidal 
mixtures. I25L l2(jt M oreover, the theory agrees well, without fitting parameters, with data from UA-MD computer 
simulations |l3l ll^j j for polyolefin blends with different local stiffness, architecture, and blend composition. 

Good agreement between the theory of mesoscopic colloidal particle mixtures and microscopic UA computer simu- 
lations of polymer blends supports the validity of our analytical renormalization procedure for polymer mixtures. The 
theory makes predictions on the evolution of the effective interpolymer potentials and chain packing as a function of 
temperature and blend composition, as well as polymer local structure, semiflexibility, and degree of polymerization. 

This paper is divided in the following way: in Sec. II, we derive monomer pair distribution functions for a binary 
asymmetric polymer blend where differences in bond length, degree of polymerization, and composition are explicitly 
taken into account. This is followed by a comparison with computer simulation data. In Sec. Ill, we present a 
renormalized description of the polymer fluid by deriving c.o.m. pair distribution functions. A brief discussion of the 
associated mean-force potential at chain-overlap is given in Sec. IV. We then show how our formalism can readily 
be cast into the language of soft-colloidal systems in Sec. V. In Sec. VI we summarize model calculations for our 
renormalized description of a binary blend and make comparisons with simulation data. A brief discussion concludes 
the paper. 



II. MONOMER PAIR DISTRIBUTION FUNCTIONS FOR ASYMMETRIC POLYMER BLENDS 



In this section we develop the theory for monomer pair distribution functions in an asymmetric polymer blend, which 
will be inputs to our renormalization approach. The theory is then tested against data from UA-MD simulations. We 
include the effect of microscopic and macroscopic parameters through a generalized form of the monomer- monomer 
pair distribution function, as derived from the thread limit in the Polymer Reference Interaction Site Model (PRISM) 
liquid-state theory by Curro, Schweizer, and coworkers . [2H Esl The choice of the theory representing the monomer 
pair distribution function is arbitrary, and other analytical approaches could have been chosen as well.@, The 
original thread model is implemented here to include the dependence on asymmetries in semiflexibility following a 
somewhat different procedure than in the original work.[30j No closure approximations are adopted in this derivation, 
and the miscibility parameter \ enters as an input to the theory, defining the proximity of the system to its demixing 
transition (the spinodal temperature). In principle, this approach can describe systems with either Lower (LCST) or 
Upper (UCST) Critical Solution Temperature. 

In the thread model, a single chain is represented as an infinite thread of vanishing thickness with hard core monomer 
diameter d — > and segment number density p — > 00, while pd 3 is finite and of order unity. Although the number 
of segments is infinite and their size approaches zero, the overall polymer dimension, given by its radius-of-gyration 
R g is finite, and enables a direct comparison with data from computer simulations of finite-size polymer chains. This 
simple chain model gives a very rough description of the liquid on the local segment scale, since it completely averages 
out solvation shells in the pair distribution function. However, this approximated description of the local structure 
becomes correct in the limit of long polymer chains, where the fine structure of the distribution function does not 
appear. 

The thread model is the analog of the model investigated in field theory approaches. Although analytically tractable, 
it has proven to be quite successful in describing properties of polymer melts, block copolymers, and blends . pll l3ll l32] | 
The main reason for its success is that it correctly captures the onset of "correlation hole" effects[33| in the monomer 
pair distribution functions at the lengthscale of the overall polymer size R g . The presence of the "correlation hole" is 
a trademark of polymer fluid structure and governs the physics of the system at the chain lengthscale. Since in our 
mesoscopic description of polymer liquids R g is the lengthscale of interest, the thread model appears to be adequate 
for the purpose of the theory developed here. 

We start from a binary blend of A and B homopolymers, having degrees of polymerization Na and Ng, and unit 
lengths oa and <tb, respectively. The polymer volume fraction is cj) — haNa/ '(uaNa + ubNb), where ua is the 
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number of molecules of type A in the mixture. Non-bonded interactions are quantified by the single \ parameter 
which describes the "monomer" interchange energy for the blend. The generalized Ornstein-Zernike (OZ) matrix 
relation for total site-site correlation functions of a binary fluid mixture in Fourier space[]3 Ei3 is 

HO) =n(jfe)C(fc) [tl(k) +H(fc)] , (1) 

where the matrices in Eq.Q are block matrices of rank 2. Here, H a p(k), C a ^(fc), and fl a p(k) are Fourier transforms 
of the corresponding correlation functions; for example, 

H afj {k) = — rH a0 (r)sm(kr)dr . (2) 
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In real space, H a p{r) = p a pph a p(r) is the (chain-averaged) site-site total distribution function matrix; C a p(r) is the 
intcrmolccular direct correlation function matrix; and O a p(r) — p a u> a (r)S a p is the intrachain structure factor matrix, 
with a, j3 € {A, B} and p a = n a N a /V the number density of monomers a inside volume V. The site averaged 
intrachain correlations are described by the Pade approximant to the Debye function (which introduces only a 15% 
maximum error to the exact expression |3^) 

Although approximated, Eq.Q enables an analytically tractable solution to the intermolecular c.o.m.c.o.m. potential, 
which is the ultimate goal of the theory developed. Moreover, the use of the Pade approximant gives good agreement 
with simulations for the total pair distribution function in both real and reciprocal spaces(2^. Agreement between 
the Debye and the Pade forms improves as N increases. 

Partial structure factors are derived from Eq.Q and S(k) = fl(k) + H(fe) as[27j 

S AA (k) = p4>u A (k)[l - p{\ - 0)w B (fc)]/A(*O 

S BB {k) =p{l-4>)u B {k)[l- pct>u A {k)]/k{k) (4) 
Sab (A) = P 2 0(1 - (j>)w A [k)u] B {k)C AB {k)/k{k) 

with p the total segment number density of the blend and 

A(fc) = l - p<jxj A (k)C AA (k) - p(l-4>)u B (k)C BB (k) + 

+ p 2 4>{\ - 4>)iv A (k)uj B (k)[C AA (k)C BB (k) - C\ B {k)\ . (5) 

The direct correlation function C aj a(fc) represents the fluid-averaged intermolecular pair potential, which is short- 
ranged or independent of k for ka -C 1. Since the lengthscale of interest for a mesoscopic description (R g ) is larger 
than er, the direct correlation function can be assumed to be independent of k, 

C aP (k) « C Q/3 (fc = 0) = Cf -e a0 /p. (6) 

The first term in Eq.©, Cg < 0, is the hard core repulsion of particles at direct contact (the potential (fcsT)~ 1 t;(r) as 
— C(r)), while the second term, e a p/p < 0, defines the "tail" interaction potential between species a and f3. 
governs the pure athermal packing of monomers in the blend, so it contains the entropic contribution that can lead 
to demixing in the high-temperature region. Consistently with the thread model we assume that the hard core 
components of the direct correlation function become as Co for any combination of A and B. Introducing Eq.Q 
and Eq.@ in Eqs.Q enables the factorization of partial structure factors into two separate contributions related to 
density and concentration fluctuations 



S AA {k) 



12p(/) 2 [\-4> 1 7 



2 



AB 



fc 2 + a 2 k 2 + c P 



-2 



S AB (k) 



12^(1-0) f 1 1 \ 



AB 



k 2 + C 2 k^+Cp 



4 



where a 2 AB 



t>a g + (1 — 4>)cr 2 A . Density fluctuations are characterized by the correlation length 
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Concentration fluctuations follow an incompressible-like relation Saa = Sbb 
the correlation length 
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-Sab, |33j and are characterized by 
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which diverges at the spinodal temperature. Here, x/p is the analog of the Flory-Huggins interaction parameter given 
by X = e AB — (caa + €bb)/2. At the spinodal temperature \ ~ * Xs, where Xs — [2Na4>] _1 + [2iVs(l — The 
ratio between A and -B chain lengths is given by p = Na/Nb, while 7 = <tb/cta defines the ratio of local stiffness 
between species A and B. Since the definition of 7 is arbitrary, for convenience we assume that the stiffest component 
in the mixture is B. Also, 
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is the correlation hole lengthscale for polymer a. Eq.© implies that asymmetry in the segment length between two 
components (7 ^ 1) increases the range of concentration fluctuations. 

The total site-site correlation function is derived from the partial structure factors as 



p 2 <P 2 h^(k) 
p 2 (l~0) 2 h^(k) 



SAA(k) - p<f)U A (k) 

Sbb (k) - p(l - <j>)uj B {k) 
P 2 4>{l-4>)h^{k) = Sab {k) , 



(11) 



where we introduced the label mm to distinguish monomer correlation functions from those of the c.o.m., which 
are derived in Section III. Eas. l|ll|l reduce in real space to the generalized thread expressions for the monomer pair 
distribution functions of asymmetric polymer blends 
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Eas. l|12|l describe the distribution of elementary units in a binary blend as a function of the intramolecular structure 
and thermodynamical parameters. In these expressions, the three lengthscales of interest £ p , £ c , and £^ appear 
uncoupled. Even in the athermal regime where x ^ Xs the mismatch in size (p. ^ 1) and/or in local semiflexibility 
(7 7^ 1) between the two components can lead to non-uniform mixing, corresponding to /i^a(^) 7^ ^ss( r ) hAB{r)- 
However, if the system is totally symmetric (Na = Nb and a a = cb), Eas. (|12fl recover the PRISM equations for a 
totally symmetric blend, which for athermal conditions recover the melt thread- model equation. |27| 

The equations just derived cannot be compared directly with simulation data of monomeric pair distribution 
functions since the thread model does not describe the presence of solvation shells in g(r) = h(r) + 1 as observed in 
simulations. However, a good estimate for the accuracy of the derived expressions can be achieved by comparing with 
simulations the number of particles included within a sphere of radius r, as 

(13) 



np(r) - S afj = 4ttp/3 / r g al 3(r)dr 
Jo 

with p/3 £ {(f>p, (1 — <p)p} for (3 g \A } B} , respectively. Tests of our monomer pair distribution functions are data from 
UA-MD computer simulations. [l3l fl4f The simulation procedure has been described in detail in recent papers^J 
Il4j and will not be discussed here. Systems investigated are blends of isotactic polypropylene (iPP), head-to-head 
polypropylene (hhPP), polyisobutylene (PIB), syndiotactic polypropylene (sPP), and polyethylene (PE). Simulations 
and blend parameters are described in the Table. 

In the simplified theoretical framework of the previous section, we enforced the limit of large polymer chains, 
which leads to a single value for the density fluctuation lengthscale, independent of the blend component. Here we 
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re-introduce finite-size effects, local semiflexibility, and branching, which are specific to each component, through 
a melt-like description of the local density fluctuation lcngthscale as £_~^ — npa^/3 + and a 6 {A 7 B}. For 

the cross pair distribution, we have £Z~g = npa AB /3 + The statistical segment length is calculated from the 

polymer radius of gyration as a a = R g A y 6/-/Va and <tb = R g B yW^s- These formulas hold in the limit of long 
flexible chains obeying Gaussian statistics. Systems investigated here include linear and branched polyolefins with 
small, densely-packed pendant groups. It is known that, in general, hi ghly branched polyolefins are more flexible 
than their linear counterparts having the same degree of polymerization J^Zl Moreover, blends of polymers with 
closely-spaced and/or small-sized pendan t group s are characterized by less efficient local packing of units. These 
effects produce a locally-disordered liquid, 39, 40] consistent with the "random packing" of a Gaussian chain and with 
the smooth shape of the thread monomer distribution function adopted here. 

Most blends considered in the simulations are far from their critical temperature. For some samples, an equation for 
the x parameter is reported in the literature, as summarized in the Table, and has been used here for the calculation 
of the concentration fluctuation correlation length. However, the size of effective sites considered in the analysis of 
the experimental data is in general different for the two components, a a &b- In these cases we normalize the value 
of x by the average number of sites per monomer |l4j| so as to be consistent with the Flory-Huggins spinodal equation 
adopted here. 

To describe samples for which an equation of x is unknown, we perform calculations in the athermal limit (x ^ Xs) 
where concentration fluctuations are minimized. Most of the blends considered here have cross contributions AB 
which could be approximated fairly well by an arithmetic or a geometric average of the self terms, AA and BB. 
This is a characteristic property of homogeneous liquids, and for blends it occurs in the athermal regime where the 
components are randomly mixed (if entropic packing effects are not dominant). 

Using the procedure just described, we calculated the number of units included in a sphere of radius r around a 
tagged polymer a. The theory shows good agreement with simulation data already at an intermolecular distance of 
a few site lengths (Fig.l). While the thread model cannot reproduce simulation data for r — > 0, it correctly describes 
the average liquid structure at intermediate and large (R g ) lengthscales, which are the lengthscales of interest in our 
mesoscopic description. We observe that the specific PRISM thread model adopted here tends to overestimate the 
number of self and cross contacts. For a UCST blend, the self-correlation functions AA and BB are best described by 
the equations in athermal conditions (chi << chi s ), while the cross-correlation contribution AB is best represented by 
the equation in thermal conditions. These trends are visible also in c.o.m. total pair distribution functions presented 
in Section VI, a feature that appears to depend on the choice of monomer pair distribution functions. 



III. CENTER-OF-MASS SOFT-CORE POTENTIAL IN POLYMER BLENDS 



To derive an analytical expression for the c.o.m. intermolecular distribution functions in a polymeric mixture, we 
extend to polymer blends a procedure outlined by Krakoviack, Hansen, and Louis |4lj for homopolymer solutions. 
While the contribution due to real sites (monomers or effective units) is averaged in the usual PRISM-like fashion, 
the c.o.m. is included in the OZ relation given by Eq.Q as an effective "auxiliary" site. For a binary fluid mixture, 
the matrices in Eq.QJ become block matrices of rank 4. We define the site-site total distribution function matrix as 
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In real space H™™^) = p a pph', 
p aPch ,ph™C{r) & ndH^{r) 
density of chains of species a; finally, H^'F(r) 
The direct correlation function matrix, 



(r) is the (chain-averaged) site-site total distribution function matrix; H£jf (r) 



Pch.aPph'^'T'ir) are the site-c.o.m. total distribution functions where p c h,a is the number 
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Pch,aPch,/3h a f}{r) is the c.o.m. -c.o.m. total distribution function matrix. 



C(k) = 



C mm (k) 




(15) 



includes the condition that direct correlation functions between auxilary and real sites as well as between two auxiliary 
sites are negligible, while the "monomer" direct correlation functions are defined in Eq.©. The intramolecular 
distribution function matrix is given by 



n mm (k) 
n Cm (k) 



n mC {k) 

Pch 



(16) 
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with fl^(r) 



mm 



[ r )S a the site-site intrachain structure factor matrix, and (r) 



Pch, a u™ c \r)S af 3 the 



site-c.o.m. intrachain structure factor matrix. The block in the upper-left quadrant in each matrix defines correlations 
between "monomer" units, whereas the block in the lower-right quadrant defines correlations between centers of mass. 
The two remaining off-diagonal quadrants contain information on the correlation between real and auxiliary sites. 
By using the definition of the static structure factor S(k) = fl(k) + H(fc), the "monomer" structure factor matrix 
recovers Eas. ()ll|) as it should. It also correctly reproduces Eqs.JHJ when the relation S(k) = [1 — D,(k)C(k)]^ 1 fl(k) 
is enforced. If the two species in the binary blend are assumed to be identical, Eas. (|14ll6[) recover equations for 
a polymer melt. [53 In general we have that h^jp(k) = h^(k) while h%™(k) ^ h™^(k) with a/(3. Using these 
conditions, Eqs.Q are solved to obtain the c.o.m. total correlation functions in Fourier space 



u™ c {k)u^ c (k) 
w2 im (fc)w™ m (fc) 
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Ea. (|17Jl formally connects c.o.m. distribution functions to the "monomer" intra- and intermolecular distribution func- 
tions, so that it is possible to calculate mesoscale properties from monomeric-scale information. For the intramolecular 
structure factor u>™ m (k), we adopt the Gaussian intrachain distribution of Eq.Q. The "monomer" -c.o.m. intramolec- 
ular structure factor (fc), can be approximated well in fc-space by a Gaussian distribution|42j as 



^ c (k)=N a e- k2R U\ 



(18) 



with a S {A, B}. Including Eas. HH3ll8fl into Eas. (|17|) yields the following expressions for the total correlation 
functions between the c.o.m. of two chains in the mixture 
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where I^(t) and I^pir) identify the concentration and density fluctuation contributions, respectively. We introduce 
here a compact notation with the function I^g(r) defined as 
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where £\ G £p} an d C P — I po\ B ) . Radii of gyration in the blend are defined as R ga f} = \J {R ga + i?^g)/2 = 

^cq/3v2 so that, for example, if a — (3 — A then R 9 aa = R g A = £, c aV% as in Eg. lfTUI) . 

Eas. p9(l can be further simplified in particular cases. For example, the density fluctuation contribution in both self 
terms is formally identical to the total distribution function between the c.o.m. of two interacting polymers in a melt 
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This expression can by reduced to a more compact form in the limit of large molecules (N — > oo) by expanding the 
total correlation function as a function of the vanishing parameter £ p / ' R g — + 
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(24) 



Since our initial equations of intra- and intermolccular structure factors rely on the assumption of Gaussian-chain 
statistics, a condition which formally holds only for polymer chains of infinite length, Eq. l|24[) is entirely consistent 
with the general description of our system and can be adopted in general for N > 30. 

The density fluctuation contribution due to cross interactions AB cannot be simplified to the melt expression, 
unless the blend is totally symmetric (a a = o~b and Na = Ng) or only slightly asymmetric. In the totally symmetric 
case, I AA = l p BB = I AB = h(r). If the system is only slightly asymmetric, the condition "Oabi ~ $AB2 is still fulfilled, 
but only for large chains (N > 30) and fi = 0.1 — 2 while 1 < 7 < 2. For this case, Eas. (|19|l simplify to 



$ a pxI p {r,£ p /R ga p^Q>) 



(25) 



with IP(r,t p /R gaf 3 -> 0) defined by Eq.©. 

For the concentration fluctuation term, the problem is further complicated by the fact that the functions in Eas. i|19|) 
depend not only on the parameter £, P /R g — > 0, but also on the additional parameter ^/R g which does not approach 
zero either in the athermal (£^ — > £ c ss R g ) or in the spinodal (£^ — > 00) regions. However, at the spinodal, where the 
concentration fluctuation lengthscale diverges, the expansion of I^pi 7 ") i n the limit of R ga p/£</> —> yields: 
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where we used the fact that lim^-^ d a pi = 1 with i E {1, 2}. The first two terms on the r.h.s. of Ea. (|26|l describe 
the spinodal decomposition and correspond to the maximum possible contribution to the pair distribution functions 
due to concentration fluctuations. When the system approaches its spinodal decomposition, contributions in R ga p/£,(t> 
become increasingly small in magnitude, to the point of no longer compensating the first two contributions on the 
r.h.s. of the equation. In the self terms of the total correlation function, this phenomenon manifests itself as an 
overshoot in the function, which corresponds to clustering of like species, or attractive interactions between molecules 
of the same type. Meanwhile, the cross correlation decreases, showing the beginning of phase separation (refer to 
model calculations in Sec. VI). 

In the athermal limit where x/Xs — * 0, if the blend is asymmetric, entropic contributions can induce demixing 
for the mesoscale liquid structure both through density and concentration fluctuation contributions. However, if the 
blend is symmetric (a a = o~b and Na = Nb), only density fluctuation contributions (I p ) govern the properties of the 
mixture, and I AA = I BB = lt R — 0. For symmetric blends, the total correlation functions recover the expression for 
the melt correlation function. [23j 



IV. CHAIN-OVERLAP AND INTERMOLECULAR MEAN-FORCE POTENTIAL 

The overlap value of intermolecular total pair distribution functions provides information on the degree of inter- 
penetration between two chains as a function of monomer density, temperature, blend composition, chain degree of 
polymerization, and asymmetry ratios (7 and /1) . In general, the intermolecular potential of mean force for two parti- 
cles in a liquid is given by W a p(r) — —kBT\ng a p(r) w —kBTh a p(r). For polymer melts, the (repulsive) potential at 
contact WTO) is known to decrease with increasing degree of polymerization, increasing density, or increasing polymer 
stiffness. |l8l l20j For polymer blends, the situation is more complex since the mismatch in local structure and chain 
length between the two components can give relevant concentration fluctuation contributions even in the athermal 
regime. 

At contact, the blend density fluctuation term is similar to the melt case. In the limit r — * and N — > 00, Ea. l23|) 
reduces to 

39 F\ 3 

lirn ^(r) « J = I p (0^ p /R gaP - 0) (27) 

■r^o 16 V 7r irpo- AB R ga p 

with i!} a pi = 1 and i £ {1,2}. Ea. (|27|) is valid for self correlation functions, and for the cross correlation functions if 
the system is symmetric or slightly asymmetric. The concentration fluctuation term in the limit r — > also simplifies 
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to 



lim I a {r) 



3 3 3 

-7\ -V a /31 2 5 

4 V 7T np(T\ B R gal 3 



e 



cot {3 



ca.0 

IT 



R, 



( 1 R 



ga/3 



(28) 



Summarizing, the total correlation functions at contact for structurally symmetric ($abi = 1) arid slightly asymmetric 
($abi ~ 1) polymer blends are written in short-hand notation as 



lim h AA {r) 

lim h B B{r) ~ 
lim hAB{r) 



1 



-7^(0)+ 7^(0,^/^-0) 



1 - 



;t(0) + f 2 / p (0,WV-0) 



(29) 



I%(0) + $ABiI p (0,Z P /R gAB - 0) 



For both density and concentration fluctuation contributions, each total distribution function at contact depends on 
density, stiffness, and degree of polymerization as I p '^ oc l/(pa\ B R ga p), essentially recovering the melt behavior: an 
increase in polymer-polymer contact is observed with increasing degree of polymerization (i? g ), stiffness (cr and R g ), 
and liquid density (p). 



V. MAPPING BLEND CORRELATION FUNCTIONS ONTO A MESOSCOPIC SOFT-COLLOIDAL 

DESCRIPTION 

Our analytical c.o.m. correlation functions effectively map a polymer blend onto a fluid of interacting soft-colloidal 
particles. To test the quality of our renormalized description we compare its predictions with known properties of 
colloidal mixtures. |34| 

At the mesoscale, a polymer blend is a two-component mixture of colloidal particles of type A and B with volume 
fraction 4> = ha/ in a + ns), total density p c h — p/N, and particle size R g A and R 9 b- When species A is chosen to be 
our reference system, the size mismatch parameter is 7 = R g B/RgA, and the reduced chain density is p* h — p c h,R? g A- 
This renormalized description can be formally obtained from Eas. df 7fl by setting the effective number of statistical 
segments in the two components to be equal, Na = Nb = N, while the chain asymmetry is completely accounted for 
by the different statistical segment lengths a' A = R g A\/&/N and a' B = R g B\/Q/N. Each effective segment includes 
the effect of branching and chain semiflexibility, and must be equal or larger than its polymer persistence length, a 
condition fulfilled by long semiflexible chains. 

It is convenient to define universal renormalized quantities which simplify the final equations. The normalized 
density fluctuation correlation length is £ p — t; p /R g A, the normalized concentration fluctuation correlation length is 

U =t*/R B A = {[1 - 0(1 - 7 2 )]/[2(l - X/Xs)}} 1/2 , and £' p = Z' p /R gA = {2^^ + (1 - tffr 2 )]}" 1 . The total pair 
distribution functions Eas. l20f) now depend only on mesoscopic reduced variables, including the space coordinate 
f = r/ R g A as 

jA (f) = 1 ML^^fW).!! fi-iVke-Vp^ 
Q/3W 4 V 7T a \ 2eJ 2 f I 2£2 1 " 



(i r -n/3 \ _ f it I a f\/3 



f /^erfc -=—= + — - e- f /^erfc 



ixV3 2a J U A V3 2a 



(30) 



where £a € {£ P , Here a = l if a = /3 = A, a = j if a = fl = B, and a = y/ (1 + 7 2 )/2 for the cross terms. Also 
&l = 62 = 1 for the self terms, while for the cross terms 

bl ~ el-d + 7 2 )/4 ' (31) 



and 

_ (2a-l)(2^- 7 2 ) 
2 [2C 2 -(l + 7 2 )/2] 2 



(32) 
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To make contact with the theory of colloidal particle mixtures (e.g., liquid alloys) it is convenient to analyze the 
properties of the renormalized polymer blend in reciprocal space. The static structure factors for the mixture are 
defined as 

SAA(k) = <p + <p 2 p ch h AA (k) , 

S BB (k) = 1 - $ + {I ~ 4>? Pchh BB {k) , (33) 
SA B (k) = 0(1- 4>)pchhAB{k) . 

Linear combinations of these functions describe fluctuations in number density and concentration, following Bhatia 
and Thornton's formalism. |25j The density fluctuation contribution (S NN in the conventional notation for "metal 
alloys") is given by 

S pp (k) = S AA (k) + S BB (k) + 2S AB (k) , (34) 
while the concentration fluctuation contribution (S ) is 

S++(k) = (1 - cb) 2 S AA (k) + <p 2 S BB (k) - 20(1 - cb)S AB (k) , (35) 
and their coupling (S CN ) is 

S p *{k) = (1 - 4>)S AA (k) - cbS BB (k) + (1 - 2cf>)S AB (k) . (36) 

Analytical formulas for Eas. i|34l3tj|) are readily obtained using Eas.l|17|). Their model calculations are presented in 
Fig. 2. Eas. l|34l36|) follow closely the behavior in reciprocal space observed for colloidal mixtures. For example, 
S pp {k) has a fc-dependence similar to the static structure factor for a single-component liquid. However, since our 
colloids are soft and can interpenetrate, there is no formation of solvation shells in the mixture and S pp {k) becomes a 
monotonically-increasing function of k. S^(k) and S p ^(k) oscillate about the values 0(1 — 0) and zero, respectively, 
as observed in colloidal mixtures with oscillations becoming less pronounced in S^(k). A more intuitive picture of 
the density-concentration fluctuation coupling term can be established by rewriting it as|25j| 

S*'(k) = 0(1 - 4>) Poh J [P A (r) - P B {r)f^^A^r 2 dr (37) 

where P a {f) = (1 — <j))gBa(f) + <j)g Aa {r) is the probability of encoutering particle clustering of species A or B around 
the colloid a E {A,B}. In this way, the function S^ p (k) represents a measure of the difference in local clustering 
between species A and B. Maxima and minima in the function point at lengthscales characterized by large asymmetry 
in the liquid structure due to the mismatch in particle size. If the two species are identical, the mixture is uniform 
and S^ p {k) = for any k. 

In the k — * limit, the density fluctuation contribution and its coupling with concentration fluctuations reduce to 
the simplified forms 

c2 



^ 07 2 + l 

2 

cA 



S PP W = jf ^' V " , (38) 
&A T 

5^(0) = 0(1- 0)2^_1 JfL , (39) 



^cA 



while the concentration fluctuation contribution is 



2 (i-0) 2 ( 7 2 -i) 2 e 



2 



S^(0) = p -fL + Y . Y \ K ' / . (40) 

1 - x/Xs (07 2 + 1 - 0)7 2 Ha 

For a blend of symmetric polymers, where R gA = R gB and 7 = 1, our equations become completely consistent with 
the theory for a mixture of symmetric colloidal particles: [25l l34l| Eci. (|38J) simplifies to S pp (0) — (^ p /^ c a) 2 the melt 
compressibility, S fp ^(0) = 0, and the concentration fluctuation contribution becomes 

s *» - .-ffr-tU ■ (41) 

Here we introduced Flory's definition of the spinodal Xs, and the renormalized x parameter for the coarse-grained 
polymer mixture AE — Nx = Ne AB — (Ne AA + Ne BB )/2. For long polymer chains or high density, Eas. l|38l39|) 
vanish since £ p /£, c — ► 0, while the concentration fluctuation contribution yields 5^(0) = 0(1 — 0)/(l — x/Xs)- 
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In general, asymmetry between the two colloidal species is quantified by the dilation fact or [2 1 



S = 



<t>v A + (1 - <t>)v B 5*^(0) 
(0 7 2 + l-0)(7 2 -l) 

Mi* -VP + 



(42) 



If the partial molar volumes per particle v a = (dV/dn a ) n g^ at p t T of the two species are identical (6 = and 7 = 1) there 
is no correlation between the fluctuations in particle number and concentration and S p ^(k) = 0. At the spionodal, 
S^(k) diverges and S -► 0. 

The isothermal compressibility for a colloidal mixture is defined [2(| as 



Pchk B TK T = S pp (0) 



S p ^(0f 



Here, 



PchkBTKT 



£ <h 2 + 1 



^ 2 



cA 



r 



(7 2 - 1) 2 ^(1 



(l-0)(7 2 



1) 2 /7 2 +^p 



(43) 



(44) 



recovers the melt compressibility when the system is composed of colloidal particles of identical size, 7 = 1. The 
compressibility is slightly temperature-dependent through the correction contribution due to S p ^(0) 2 / S^(0). The 
latter, however, is small for large polymer chains since it scales with degree of polymerization as iV -1 in the athermal 
regime, and vanishes approaching the spinodal curve where the concentration correlation length diverges. This is true 
at any lengthscale and the ^-dependent blend compressibility can be approximated by its first contribution S pp (k) for 
the entire range of k and for blends of polymer chains with symmetric or asymmetric size, in agreement with colloidal 
particle mixtures (Fig. 2). 



VI. CENTER-OF-MASS DISTRIBUTION FUNCTIONS: MODEL CALCULATIONS AND 
COMPARISON WITH UA-MD COMPUTER SIMULATIONS 

In this section we present some predictions of our theoretical approach for the clustering and interpenetration of 
polymer chains in a blend as a function of chain length, semiflexibility, polymer volume fraction, and x dependence 
by model calculations. In our model calculations, unless otherwise specified, Na = 500, p m = 0.03 A -3 , and a a = 3.0 
A. These parameters have been chosen to be consistent with the UA-MD simulations of polyolefin blends presented 
below. The chain length, however, has been chosen to be quite large to avoid finite-size effects which could veil the 
general trend of blend properties. Calculations performed "close to the spinodal condition" are for a \ parameter 
which deviates 0.01% from Xs- 

As a first study, we look at a structurally symmetric blend for which 7 = 1 and Na = Nb , while varying the blend 
composition <f> 6 {0.25,0.50,0.75}. In the athermal limit x ^ Xsi au distribution functions become identical for any 
blend composition <f>. However, when concentration fluctuations develop (Fig. 3), we observe an enhanced clustering 
and increased number of self contacts in the minority (low <fi) species, which corresponds to a strong attractive 
interaction in the mean-force potential. The majority (high <fi) species increases the number of self-contacts, but still 
resembles in shape its athermal limit. This effect corresponds to the initial stage of droplet formation in the demixing 
of compositionally asymmetric mixtures. In general, the extent of demixing increases with increasing asymmetry in 
blend composition. 

To investigate the effect of mismatch in chain flexibility, we analyze blends of polymers with identical degree of 
polymerization (Na = Nb) where component A is flexible while the stiffness of component B is increasing (increasing 
R g B and 7 £ {1.0,1.2,1.5}). The data from UA-MD simulations analyzed subsequently in this paper (see Table) 
belong to this group. In general, we observe that an increase in 7 enhances the number of BB and AB contacts and 
reduces the number of AA contacts (Fig. 4). This effect is observed in both athermal and thermal mixtures. A stiff 
chain has a larger R g than a flexible molecule with the same degree of polymerization. In this case, the number of 
chains that can occupy the volume spanned by the stiff component increases with increasing R g . Furthermore, stiff 
chains pack better at short distances due to geometric effects (stretched configurations enhance BB contacts), while 
flexible chains, which have a coiled configuration, tend to pack more efficiently at distances of order R cj a- Overall, 
these two effects lead to an enhanced miscibility in mixtures of polymers with different flexibility and same degree of 
polymerization. 
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The potential of mean-force for chains of e qual length is found to decrease with increasing degree of polymerization 
in agreement with computer simulations. |l6l Il8| The effect of mismatch in chain length for polymer species having 
the same semiflexibility is studied using a blend of polymers with Na = 500 fixed and Nb G {250,500,1000}, 
corresponding to /i G {0.5, 1.0, 2.0} (Fig. 5). Analogous effects are observed in both the athermal and the thermal limits. 
In the athermal limit, the mismatch in chain size prevents the random mixing of the two components and "entropic" 
effects appear. Trivially, as the R g of the B component increases the number of chains that can interpenetrate in 
the volume spanned by B increases. This leads to a larger number of intermolecular contacts for all pair distribution 
functions AA, BB 7 and AB. This simple trend is combined with a more subtle behavior: a mismatch in chain length 
favors interpenetration between different species, so that the number of AB contacts is always higher than the number 
of AA and BB contacts except for chains of the same length. This property becomes more pronounced when the 
difference in chain length increases. Also, short flexible chains tend to pack efficiently on the scale of their R g , while 
these effects are averaged out in large flexible chains, where "solvation shells" do not appear (a behavior analogous to 
the monomer distribution functions). Our theoretical predictions are in qualitative agreement with off-lattice Monte 
Carlo simulations of polymer solutions. [Ifij 

To investigate in more detail these effects, wc analyze polymer blends which have constant R g A and R 9 b- The 
change in chain length in the B component is correlated with a change in semiflexibility. For example, the A component 
is mixed with the component B, which is both shorter and stiffer than A (fi = 2.00 and 7 = 1.88), or with a chain B 
larger than A but with comparable flexibility (fi = 0.66 and 7 = 1.09). In analogy with our previous observation, we 
see (Fig. 6) that stiff chains tend to pack more efficiently than flexible ones at short distances, increasing in this way 
the number of BB and AB contacts. As a consequence, the number of AA contacts in a mixture of stiff B chains 
is decreased with respect to the mixture of flexible polymers. Close to the spinodal, flexible chains (^4 and B) are 
in coiled configurations which tend to pack more efficiently at the distance of their R g . Stiff chains, instead, tend to 
interpenetrate more efficiently and the number of BB contacts strongly increases in the mixture containing short and 
stiff B chains. 

Next, we compare our analytical expressions against data from UA-MD simulations of polymer blends described in 
the Table. |13Lll4| For each blend considered, the theory agrees with simulation data fairly well. However, we present 
here only few systems which are representative of our calculations. Input parameters to the theory and the procedure 
used in the calculations were discussed in Section II. 

The c.o.m. total correlation function provides an estimate of the number of molecules interpenetrating at some 
relative distance r. In all the plots (Fig. 7) we observe that chains of the stiffest component (£?), which have extended 
configurations, tend to pack at short distances r < R g s more efficiently than flexible ones. Flexible molecules, which 
have coiled configurations, show a high (low) number of "intramolecular" (intermolecular) contacts and pack most 
efficiently at distances on the order of the overall polymer size, r ps \.bR g A- The extent of intermolecular chain 
packing upon blending also depends on polymer flexibility. The theory predicts that the stiff (flexible) component 
packs better (worse), and the number of self contacts increases (decreases) when mixed with a more flexible (stiffer) 
polymer, in agreement with simulations. |12| 

In general, the agreement between theory and simulations is good with the exception of the PIB /hhPP and PIB /iPP 
blends for which agreement is only qualitative. For these systems the theory overestimates the number of intermolec- 
ular contacts. However, it is well known that PIB blends are usually immiscible blends at these temperatures and 
for these chain lengths. 0, ^| PIB exhibits very efficient intra- and intermolecular packing which leads to a thermal 
expansion coefficient and an isothermal compressibility smaller than in other polyolefin blends. |43j This behavior is 
due to strong attractive interactions between the methyl (— CH 3 ) groups, which in PIB molecules are in very large 
number (about 50% of the total number of united atoms). Q However, even for PIB blends, the theory shows very 
good agreement with simulation data for r « R g . 

Finally, in Fig. 8 we show how the theory can predict the correct qualitative behavior for blends following a LCST 
phase diagram. At high temperature the blend demixes and shows a distribution function for the AB component that 
is consistently lower than the AA and the BB contributions. At low temperature, instead, the mixing of the A and 
B species is enhanced with the AB function being always higher than its A A and BB counterparts. 

VII. CONCLUSIONS 

In this paper we present the first analytical solution for the center-of-mass total distribution functions for polymer 
mixtures. We start from first-principles liquid state theory and derive a set of equations which describe c.o.m. 
distribution functions between pairs of interacting molecules in the mixture. The mesoscale description of the liquid 
so obtained depends on local chemical parameters, such as the radius of gyration and bond length. Moreover, mesoscale 
distribution functions depend on thermodynamic parameters such as density, blend composition, and proximity of 
the system to its spinodal decomposition. 
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The analytical c.o.m. total distribution functions quantitatively reproduce data from microscopic united-atom 
molecular dynamics simulations without the need of fitting parameters. In this comparison, input to the theory are 
data from the UA-MD simulations. The excellent agreement between analytical theory and simulations in both real 
and reciprocal space is an indication that our approach provides a reliable mesoscale description of polymer mixtures. 

Moreover, our formalism recovers the known analytical functions for the fluctuations in number density and con- 
centration (structure factors) for interacting soft-colloidal mixtures. The one-to-one agreement between the two 
formalisms supports the validity of our approach as a procedure to map polymer blends onto systems of interacting 
soft-colloidal particles. 
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Table: Polyolefin blends (T = 453 K, N A = N B = 96). 
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FIG. 1: Plot of n a p(r) against r for the hhPP/PE blend. Theoretical curves for athcrmal (full lines) and thermal 
(dashed lines) conditions are compared with simulation data (symbols): AA— (circles), AB— (diamonds, inset), and 
BB— terms (squares). Arrow indicates the lengthscale R g A- 

FIG. 2: Plot of colloid partial structure factors against kR g A for 7 € {1.0 (dashed lines), 1.5 (full lines)} with 
fi = 1.0 and <f> = 0.5. Left panel. Upper portion: S"^(fc) with X /Xs & {0.0, 0.2, 0.4, 0.6, 0.8} from bottom to top. 
Lower portion: S p ^(k) with x/Xs = 0.0 (top line) and x/Xs — 0.8 (bottom line). Right panel: S pp (k) for x/Xs € {0.0, 
0.2, 0.4, 0.6, 0.8} (curves are indistinguishable in the plot). Right panel also shows S pp (k) — S 2 S^(k) (dot-dashed 
lines) for 7 = 1.5. 

FIG. 3: Plots of h a p(r) against r/R g A for a structurally symmetric blend (/1 = 1.0 and 7 = 1.0) close to (0.01% 
from) its spinodal decomposition. Left panel: effect of blend composition for AA—tara at <j> e {0.25 (long-dashed 
line), 0.50 (dot-dashed line), 0.75 (full line)}. The behavior for BB— term is complementary to the AA— term at 1 — 0. 
Right panel: AB— term for which curves superimpose onto the full line. Also shown for comparison is the athermal 
case at <p — 0.5 (dashed lines). 

FIG. 4: Plots of h a p{r) against r/R g A for 7 e {1.0 (dashed lines), 1.2 (full lines), 1.5 (dot-dashed lines)} with 
p = 1.0 and 4> — 0.5 for athermal conditions. Left panel: AA— (thin lines) and BB— terms (thick lines). Right panel: 
AB— terms. Note that AA = BB for the 7=1 case. 

FIG. 5: Plots of /i Q( g(r) against r/R g A for /1 e {0.5 (dot-dashed lines), 1.0 (dashed lines) 2.0 (full lines)} with 
7 = 1.0 and (f> = 0.5 under athermal conditions. Left panel: AA— terms; middle panel: AB— terms; right panel: 
BB— terms. 

FIG. 6: Plots of h a p{r) against r/R g A for systems with different /1 and 7, but constant radii of gyration R g A 
and R g B at <j> = 0.5 for athermal conditions: full lines for p = 2.00 and 7 = 1.88; dashed lines for p = 1.00 and 
7 = 1.33; dot-dashed lines for \i = 0.66 and 7 = 1.09. Left panel: AA— terms, middle panel: AB— terms; right panel: 
BB— terms. 

FIG. 7: Plots of h a p(r) against r/R g A for blends. Theoretical predictions in athermal (full lines) and thermal 
(dashed lines) conditions are compared with UA-MD simulations of blends (symbols): AA— (circles), AB— (dia- 
monds), and BB— terms (squares). Blend parameters are given in the Table. 

FIG. 8: Plots of h a[) {r) against r/R gA for a LCST blend described by X = 0.01125 - 4.75/T (N A = 96, p m = 0.034 
A~ 3 , a a = 2.44 A, 4> = 0.50, and Xs = 0.021). Shown are AA- (dot-dashed lines), AB- (full lines), and SB-terms 
(dashed lines) at two temperatures: T = 150 K (thin lines) and T = 10,000 K (thick lines). 



